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ABSTRACT 



Context. Regarding the significant interest in both dark matter and the application of MOND to early-type galaxies, we investigate 
the MOND theory by comparing its predictions, for models of constant mass-to-light ratio, with observational data of the early-type 
galaxy NGC 4649. 

Aims. We study whether measurements for NGC 3379 and NGC 1399 are typical of early-type systems and we test the assumption 
of a Newtonian constant M/L ratio underlying most of the published models. 

Methods. We employ the globular clusters of NGC 4649 as a mass tracer The Jeans equation is calculated for both MOND and 
constant mass-to-light ratio assumptions. Spherical symmetry is assumed and the calculations are performed for both isotropic and 
anisotropic cases. 

Results. We found that both Jeans models with the assumption of a constant mass-to-light ratio and different MOND models provide 
good agreement with the observed values of the velocity dispersion. The most accurate fits of the velocity dispersion were obtained for 
the mass-to-light ratio in the fi-band, which was equal to 7, implying that there is no need for significant amounts of dark matter in the 
outer parts (beyond 3 effective radii) of this galaxy. We also found that tangential anisotropies are most likely present in NGC 4649. 
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1. Introduction 

The elliptical (early-type) galaxies are pressure-supported rather 
than rotationally supported systems and it is observationally 
more difficult to infer the presence of dark matter than in the 
more well-studied late-type galaxies. This is because, at large 
galactocentric distances, ellipticals do not have a straightforward 
tracer of circular orbits that is similar to the neutral hydrogen 
used to study the halo kinematics of spirals. Fortunately, how- 
ever we are able to use several alternative tracers in the nearby el- 
lipticals: 1) integrated stellar spectra (although the observations 
are limited practically within ~ 3 - AR,,, where is effective 
radius), 2) X-rays 3) planetary nebulae (PNe) and 4) globular 
clusters (GCs). 

Studies of ellipticals have indicated that within ~ 2 - 3Re 
dark matter does not appear to dominate (e.g. Samurovic & 
Danziger 2005, 2006). Beyond ~ 3Re, dark matter appears to 
start to play important dynamical role. Although the presence 
of dark matter in the Universe is widely accepted, some alterna- 
tive opinions have been expressed such as the theory of MOND 
(Milgrom 1983; recent reviews in Scarpa 2006, Milgrom 2008). 
The case for giant ellipticals has not been considered adequately 
and one possible reason is the scarcity of available kinemati- 
cal data out to large galactocentric distances. Several attempts to 
model ellipticals with dark matter and/or MOND theory have ap- 
peared (e.g. Schuberth et al. 2006, Klypin & Prada 2007, Tiret et 
al. 2007, Richtler et al. 2008) with ambiguous results. Notably, 
Tiret et al. (2007) claimed that a MOND model for the dynamics 
of NGC 3379 reproduced the observations on all scales, while 



Richtler et al. (2008), studying NGC 1399, reached the opposite 
conclusion - that the best-fit MOND model still requires an "ad- 
ditional hypothetical dark halo". The latter result, for matter on 
small scales, agrees with observations implying the need for dark 
matter at the centers of clusters: for example, Angus, Famaey & 
Buote (2007) found that a hidden mass component (about 1 .5 to 
4 times more massive than the total visible mass) is required to 
explain the hydrostatic equilibrium of clusters with temperatures 
ranging from 0.7 to 8.9 keV. 

In this paper, we discuss the early-type galaxy NGC 4649 
(M60) belonging to the Virgo cluster, for which new observa- 
tional data have become available (Lee et al. 2008a). This galaxy 
is more typical than either a field giant elliptical (NGC 3379) or 
the central galaxy of a cluster (NGC 1399), since it is a typical 
cluster member. The plan of the paper is as follows: in Sect. 2, 
we present the data related to NGC 4649; in Sect. 3, we calcu- 
late the total mass of this galaxy for Newtonian gravity and the 
MOND approach and solve the Jeans equation for the constant 
mass-to-light ratio and MOND approach to determine the best- 
fit parameters. Our conclusions are presented in Sect. 4. 



2. Observational data 

The data used in this paper are based on observations of GCs in 
NGC 4649 presented by Lee et al. (2008a). The sample consists 
of 121 GCs (83 blue and 38 red GCs). In all of our calculations, 
we always consider the entire sample, i.e. we do not divide the 
data into subsets, so that we are able to have a larger number of 
GCs per bin in our analyses. 
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NGC 4649 (M60) is a giant elliptical galaxy in the Virgo 
cluster that has a nearby companion, NGC 4647 (Sc galaxy at 
2. '5 from the center of NGC 4649). The systemic velocity of 
NGC 4649 is vvi = 1117 + 6 km s"'. We consider two values 
for the distance to NGC 4649: (i) the first one, d - 17.30 Mpc is 
based on the surface brightness fluctuation method and is taken 
from Lee et al. (2008a) (in this case, one arcsec corresponds to 
84 pc); and (ii) the second one is based on the aforementioned 
systemic velocity, which implies that using the WMAP estimate 
of the Hubble constant from Komatsu et al. (2008) ho - 0.70, we 
obtain d - 15.96 Mpc, for which one arcsec corresponds to 77.5 
pc. Both values of the distance are tested below and the most 
appropriate fit is presented. For the effective radius, we assume 
the value of Rg = 90 arcsec (equal to 7.56 kpc for d = 17.30 
Mpc, and 6.97 kpc for d = 15.96 Mpc), which is the value taken 
from the paper by Kim et al. (2006). We note that there are dif- 
ferent estimates in the literature, for example, Lee et al. (2008a) 
assumed Rg - IIQ arcsec and according to the RC3 catalog (de 
Vaucouleurs et al. 1991), Re = 69 arcsec. The discrepancy does 
not impact significantly our present conclusions. 

We used the radial velocities of GCs to determine the kine- 
matics of NGC 4649: we calculated the velocity dispersion, the 
skewness and kurtosis parameters, sj, and 54, using standard defi- 
nitions and the NAG routine G01AAF. We note that our intention 
is not to reconstruct the full line-of-sight velocity distribution, 
because it is known (see Merritt 1997) that for small data sets 
such as the one we consider here, that contain less than a few 
hundred objects per bin, this is impossible. We note, however, 
that Wu & Tremaine (2006) developed a maximum likelihood 
method for determination the mass distribution in spherical stel- 
lar systems that uses test particles and applied this method suc- 
cessfully to a giant elliptical M87 using 161 GCs found in this 
galaxy. 

We therefore calculate skewness and kurtosis parameters to 
determine whether, in some bins, a significant departure from 
a Gaussian distribution (and a tendency toward some particular 
type of orbits) exists. This is similar to the approach applied by 
Teodorescu et al. (2005) in their Fig. 18. Our results regarding 
departures from a Gaussian distribution are applied below where 
we apply different anisotropics when solving the Jeans equation. 
The results are given in Fig. [T] and in Table [T] What is obvious 
from the data (and what makes this galaxy a good candidate for 
the dynamical modeling based on two different techniques) is 
that the velocity dispersion appears to be constant throughout 
the entire galaxy, and, based on the values of S4 (small nega- 
tive values but consistent with zero), there appears to be a ten- 
dency towards tangential orbits; the same conclusion based on 
Jeans modeling was reached by Hwang et al. (2008) (and also 
by Bridges et al. (2006), who found that the orbital distribution 
beyond 100 arcsec becomes tangentially biased). 

For the Jeans modeling presented below, we also needed to 
determine the radial surface density profile of the clusters in 
NGC 4649. This plot is given in Fig. |2] and represents the ra- 
dial distribution of the total sample of GCs in NGC 4649 as a 
function of radius. The fit is given by a power law of the form: 
2 oc r Using the least squares method, we decided that, within 
a galactocentric radius of approximately 1 arcmin y = 0.338, 
and for the region beyond 1 arcmin, that we should use the value 
y = 1 .285 derived by Lee et al. (2008b). 
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Fig. 1. Kinematics of NGC 4649 based on the total sample of red 
and blue GCs. From top to bottom: radial velocity of the GCs in 
km s"'; velocity dispersion calculated in a given bin; and the 53 
and S4 parameters, which describe symmetric and asymmetric 
departures from the Gaussian, respectively. 
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Fig. 2. Radial distribution of GCs (both blue and red) in 
NGC 4649. A power law is fitted to the radial surface density 
of GCs: E oc r"'' (see text for details). 
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Table 1. Kinematics data for NGC 4649 for the total sample of clusters 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


232 


30 


-0.03 


0.06 


-0.04 


0.04 


30 


3 


231 


24 


0.02 


0.05 


-0.03 


0.04 


45 


5 


233 


31 


0.00 


0.07 


-0.03 


0.05 


28 


7.5 


217 


36 


-0.05 


0.08 


-0.04 


0.06 


18 



NOTES - Col. (1): central point for a given bin. Col. (2): velocity dispersion for the total sample of clusters. Col. (3): formal errors for the 
velocity dispersion of the total sample of clusters. Col. (4): parameter for the total sample of clusters. Col. (5): formal errors for the S3 
parameter for the total sample of clusters. Col. (6): S4 parameter for the total sample of clusters. Col. (7): formal errors for the S4 parameter for 
the total sample of clusters. Col. (8): number of clusters in a given bin. 



3. Models 

We model the observed velocity dispersion beyond 1 arcmin (be- 
cause our observationally based points are all beyond ~ 1 arc- 
min) using the Jeans equation (e.g. Binney & Tremaine 1987): 



dcrf. , (2/3. + a) GM{r) 

— — + cr, = 5— 

dr r 



(1) 



where cr,. is the radial stellar velocity dispersion, 
a - dlnp/dlnr is the slope of tracer density p (the sur- 
face density is given above and in the models below we use 
a - -2.285). The rotation speed v'lot was found to be non-zero, 
i.e. Vi-ot - 141^3g km s"' (Hwang et al. 2008) and is used in all 
modeling below. A parameter jS, is introduced to describe the 
non-spherical nature of the stellar velocity dispersion: 



(2) 



where Vg = vf'+crg. For < j6, < 1, the orbits are predominantly 
radial, and the line-of-sight velocity distribution is more strongly 
peaked than a Gaussian profile (positive 54 parameter), and for 
-00 < j6, < the orbits are mostly tangential, so that the pro- 
file is more flat-topped than a Gaussian (negative 54 parameter) 
(Gerhard 1993). 

In all models below, we calculated the projected line-of-sight 
velocity dispersion (e.g. Binney & Mamon 1982; Mathews & 
Brighenti 2003): 



(3) 



where the truncation radius, r,, extends beyond the observed 
kinematical point of highest galactocentric radius. We assumed 
in all the cases that r, ~ IRe- 

3.1. Newtonian mass-follows-light models 

A relation for determining the typical mass-to-light ratio in el- 
liptical galaxies was given by van der Marel (1991), who found 
for 37 bright ellipticals that the average mass-to-light ratio in the 
B-band in solar units was: M/Lb = (5.93 + 0.25) /150 becoming 
M/Lb = (1 1 .86 ± 0.50) h and therefore M/Lb = 8.30 ± 0.35 for 
h - 0.7. Since this study addressed the inner parts of ellipticals, 
we considered the absolute upper limit for the visible (stellar) 
component to be M/Lb ~ 9 - 10. An inferred mass-to-light ra- 
tio of over approximately 10 in a given region would therefore 
imply the existence of dark matter 



For constant mass-to-Ught ratio models, we consider rela- 
tions that include stellar mass and dark matter distributed in the 
form of the standard Hernquist (1990) profile: 



PH(r) = 



Mt a 



1 



27r r (r + a) 



3 ' 



(4) 



which has two parameters: the total mass Mr and scale length, 
fl, where Re = 1.8153fl. We solve the Jeans equation (Eq. [T]) 
and consider values of M/Lb of between 7 and 15, while repeat- 
ing that M/Lb^IQ is incompatible with an entirely stellar com- 
ponent. As can be seen from Fig. [3] we tested 3 values of the 
mass-to-light ratio that are presented by the hatched stripes: low 
{M/Lb = 7), intermediate (M/Lb = 10), and high (M/Lb = 15, 
includes at least ~ 50% of dark matter). In all cases, we allowed 
a realistic variation in both the parameter and the distance to 
the galaxy, d for the fixed value of the mass-to-light ratio. If we 
need to identify the constant mass-to-light ratio that most accu- 
rately describes the observed velocity dispersion data, we would 
present the model with the following parameters: lower mass- 
to-light ratio, M/Lb - 7, slightly tangential orbits, yS» - -0.2, 
and the distance d - 17.30 Mpc (see Fig. [3] and Table |2|. We 
note that a good fit can also be achieved for a higher value of the 
mass-to-light ratio but then the parameter rises and the orbits 
become isotropic (see Fig.|3] central "S2" stripe, which presents 
a case for which M/Lb - 10). 

With the assumption of the constant mass-to-light ratio the 
overall conclusion of the Jeans modeling is that, although we 
detect a hint of a rising mass-to-light ratio with increasing radius, 
due to the large error bars, we can obtain a satisfactory fit by also 
assuming a constant value M/Lb ~ 7, which implies that dark 
matter is not playing an important dynamical role, even beyond 
~ 3Re. The value of the constant mass-to-light ratio M/Lb ~ 7 
implies that the total mass of NGC 4649 is equal to ~ 5 x 10"Mo 
at ~ 3Re. This estimate is in agreement with predictions based 
on the assumption that the Virial theorem holds (see Bertin et 
al. 2002; Cappellari et al. 2006). This total mass is used below. 



3.2. MOND models 

We follow the recommendation of our referee and calculate the 
total mass of NGC 4649 in MOND gravity using three differ- 
ent formulas: (i) the "simple" MOND formula from Famaey & 
Binney (2005), (ii) the "standard" formula (Sanders & McGaugh 
2002), and (iii) the Bekenstein's "toy" model (Bekenstein 2004). 
We write the Newtonian acceleration as = aij(a/ao), where 
a is the MOND acceleration, //(x) is the MOND interpolating 
function where x = a/ao, and aq 



1.35+"-2? X lO-** cms-^ is a 



-0.42 
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universal constant (Famaey et al. 2007). The interpolation func- 
tion yu(fl/fl()) shows an asymptotic behavior, jj ~ I, for a » ao, 
and we derive the Newtonian relation in the strong field regime, 
and fi = a/flo for a <K aq. The MOND dynamical mass. Mm, 
can be expressed in terms of the Newtonian values, using 
the following expression (e.g. Angus et al. 2007): 



MM(r) = M/v(r) X fi(x). 



(5) 



The intepolation function can have different forms as given be- 
low. 



1. A "simple" MOND formula is given by: 



1 + X 



(6) 



In this case, the circular velocity curve can be written as (e.g. 
Richtler et al. 2008) 



y2 

' circM 



y2 y4 

— T — + A — : — + y X ao X r. 



where VchcN is the Newtonian circular velocity. 
2. A "standard" MOND formula is given by: 

Y 

fi(x) = 



(7) 



(8) 



It can be shown that the circular velocity curve then be- 
comes: 



v: 



circ.M 



+ V^. X at. X r^. 

circ.N 



3. Finally, for the "toy" model the MOND formula is: 



fi(x) = 



-1 + Vl + 4x 



1 + Vl +4x 



(9) 



(10) 



After some arithmetic, we derive the expression for the cir- 
cular velocity curve: 



circ.M 



^cii-cN + VaO X r VcircN- 



(11) 



These three MOND models are tested below using Jeans models, 
and the best-fit model parameters are indicated in each case. 

We solved the Jeans equation (Eq. [U based on the Hernquist 
profile (Eq. HJ, where Mj- is calculated using Eq. |5] We varied 
different parameters: distance d, parameter aq, and anisotropy 
yS*. The best-fit models are presented in Fig. |4] and Table |2] We 
tested all three MOND models and it can be seen that the best fits 
were obtained for higher value of the distance, d - 17.30 Mpc. 
We note that the quality of the fits for all MOND models was 
approximately the same, as in the case of the constant mass-to- 
light ratio models, as indicatated by the reduced values (see 
column 6 in Table|2]i. We formed our Jeans models based on an 
expression for the mass corresponding to each of the constant 
mass-to-light ratios, as in the previous section. The lower mass- 
to-light ratio {M/Lb - 7) provides a good fit to the observed 
velocity dispersion for all tested models. The preferred value of 
the constant oo - 1.35 x 10"** cm s"' was found in two cases 
(for the "simple" and "toy" model), whereas a lower value was 
found for the "standard" model, oq = 0.93 x 10"^ cm s '. In all 
tested models, we found that there is a detection of tangential 
anisotropy, the largest extent for the "toy" model (yS, - -0.5) 



and the smallest extent for the "standard" model (/?, = -0.3). 
It is important to emphasize that we cannot exclude higher (but 
not extremely high, see below) mass-to-light ratio in our models; 
for example, the "toy" model with M/Lb - 9,ao - 0.93 x 10"^ 
cm s"', c/ = 17.30 Mpc, and (J3, = -0.3) also provides a satis- 
factory fit to the observed data with reduced ;if2 - 2.90. In gen- 
eral, because of the large observational error bars at this stage, 
we cannot exclude the possibility of an increase in the velocity 
dispersion in the outer parts, which would imply the increase of 
the total mass-to-light ratio and the existence of dark matter We 
emphasize that these models are based on a constant mass-to- 
light ratio and, therefore, a fixed higher value of the mass-to- 
light ratio, in principle corresponds to a higher velocity disper- 
sion close to the center and a more realistic value in the outer 
parts ofNGC 4649. 

It is known that the "external field effect" is a phenomeno- 
logical requirement of MOND, which has strong implications 
for non-isolated systems (Sanders & McGaugh 2002), such as 
NGC 4649. This effect is complicated and was analyzed in de- 
tail by Wu et al. (2007); here, we provide only an approxi- 
mate estimate of this effect on our results. We follow Famaey 
et al. (2007b) and calculate the gravitational force per unit mass 
exerted by the Virgo cluster on NGC 4649: 



(GM^.flo) 



1/2 



^ 0.15fl() 



(12) 



where My^^. = 2.1 x 1O'"'M0 is the baryonic mass of the Virgo 
inside dc ~ 1 Mpc (which is the distance of NGC 4649 from the 
center of the Virgo cluster; the distance to NGC 4649 is equal to 
17.30 Mpc and is also used) taken from the paper by McLaughlin 
(1999) (his figure 2) who estimated that a global baryon fraction 
in Virgo is ~ 7%. We then apply the external effect by insert- 
ing /i(|a + aextl/flo) only in the case of the "toy" model to infer 
the effects on the Jeans modeling. The result is plotted in Fig.|4] 
with the dotted line (and is also given in Table |2]i. The modeled 
velocity dispersion is higher than in the case when we neglect 
the external field but not by much. This result is again consistent 
with a hypothesis of no dark matter. Parenthetically, the same 
reasoning applies much more forcefully to the acceleration due 
to the presence of the galaxy companion NGC 4647, for which 
the value in Eq. [12] is smaller for between one and two orders 
of magnitude, depending on the baryonic mass estimate for the 
companion. 

The same conclusion regarding the dynamical importance of 
dark matter reached in the case of the constant mass-to-light ratio 
also holds for the Jeans models based on MOND: successful fits 
were obtained without dark matter throughout the entire galaxy 
for all MOND approaches. 

4. Conclusions 

The GCs kinematics of NGC 4649 was studied out to 260 arc- 
sec {- 2.9Re using the value of effective radius adopted in this 
paper) by Bridges et al. (2006). They found, using isotropic and 
axisymmetric orbit-based models, that dark matter exists in the 
halo of this galaxy. 

In this paper, we have studied the dynamics of the early-type 
galaxy NGC 4649 (M60) for both Newtonian and MOND ap- 
proaches using the Jeans equation and found that there is no need 
for significant amount of dark matter in its outer parts. The prin- 
ciple results of our study were: 

1 . We studied the kinematics of NGC 4649 and found that the 
velocity dispersion remains approximately constant through- 
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Table 2. Best fits for constant M/L ratio and MOND modeling 
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Model 


MILb 


flo [10-* cm s-2] 


d [Mpc] 


p. 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Const 


7 




17.30 


-0.2 


1.36 


MOND Simple 


7 


1.35 


17.30 


-0.4 


2.31 


MOND Standard 


7 


0.93 


17.30 


-0.3 


1.66 


MOND Toy 


7 


1.35 


17.30 


-0.5 


1.83 


MOND Toy (ext. field) 


7 


1.35 


17.30 


-0.5 


1.78 



NOTES - Col. (1): the name of the model; "ext. field" is for the MOND toy model for which the external field was talsen into account (see text 
for details). Col. (2): mass-to-light ratio in the 6-band used. Col. (3): the best-fit function value of the constant a^; it is not used in the case of 
the model with the constant mass-to-light ratio. Col. (4): the best fitting value of distance to NGC 4649. Col. (5): the best fitting value of the jS, 
parameter. Col. (6): reduced of the best fit for a given model. 




Fig. 3. Jeans modeling of NGC 4649 using the constant M/L ap- 
proach. Three different stripes, which correspond to three differ- 
ent constant mass-to-light ratios, are given and the upper limit 
is always determined for yS* = -0.2, and d - 15.96 Mpc while 
the lower limit is determined for yS, = 0.2 and d = 17.30 Mpc. 
The lowest stripe "SI" is for the constant mass-to-light ratio 
M/Lb = 7, the stripe in between, "S2" is for M/Lb -10, and 
the upper stripe, "S3" is for M/Lb - 15. The solid line is the 
best fit obtained for the constant mass-to-hght ratio models: in 
this case M/Lb = 1,13* = -0.2 wAd = 17.30 Mpc. The upper 
scale corresponds to the distance d = 17.30 Mpc. 



out the entire galaxy. We detected an observational hint of 
small tangential anisotropics in this galaxy, which agrees 
with other results in the literature (Bridges et al. 2006, 
Hwang et al. 2008). 
2. We used the Jeans equation to infer the mass-to-Ught ratios 
that can describe the observed velocity dispersion profile of 
NGC 4649: we found that the successful fit to the velocity 
dispersion can be obtained using M/Lb ~ 7 with moderate 
tangential anisotropics that indicate the lack of dark matter 
throughout this galaxy. Higher values of the mass-to-light 
ratios (M/Lb ~ 9 - 10), which are still consistent with a 
no dark matter hypothesis, are also permitted but then the 
anisotropics tend to vanish. 



Fig. 4. Jeans modehng of NGC 4649 using the MOND approach 
and M/Lb = 7. Solid line is for the "standard" model for which 
flo = 0.93 X 10"** cm s"^ and y6, = -0.3. Long-dashed line is 
for the "simple" model for which ao = 1.35 X 10 cm s"^ and 
j8, = -0.4. For the "toy" model we present two fits for which 
flo = 1.35 X 10"^ cm s"^ and/3, = -0.5: (i) short-dashed line for 
the case when the external field is absent (ii) dotted line for the 
case when the external field is present. In this figure the upper 
scale corresponds to d = 17.30 Mpc because the best fits in all 
the cases were obtained using this distance. 

3. We also used the Jeans equation to test the predictions of 
the MOND theory using three different models: "simple", 
"standard", and "toy". We found that all the tested MOND 
models provided a good fit to the observational velocity dis- 
persion data (again for M/Lb ~ 7 with moderate tangential 
anisotropics). 

4. We have also briefly addressed the problem of the surround- 
ing environment, since NGC 4649 belongs to the Virgo 
cluster. We found that the external field at the position of 
NGC 4649 does not significantly alter our conclusions: the 
Bekenstein "toy" model again provides a good fit to the ob- 
served velocity dispersion (assuming M/Lb ~ 7) and is con- 
sistent with a no dark matter hypothesis. 
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